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1.  INTRODUCTION 

The  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS)  is 
the  Department  of  Defense's  (DoD's)  high- 
resolution  global  numerical  weather  prediction 
(NWP)  system.  Its  development  and  operation 
is  a  joint  activity  of  the  Naval  Research 
Laboratory  (NRL)  and  the  Navy's  Lleet 
Numerical  Meteorology  and  Oceanography 
Center  (LNMOC).  NOGAPS  is  a  complete 
operational  forecasting  system  that  includes 
data  quality  control  (Baker  1992),  tropical 
cyclone  bogusing  (Goerrs  and  Jeffries  1994), 
operational  data  assimilation  (Daley  and 
Barker  2001),  balanced  initialization  (Errico  et 
al.  1988),  and  a  global  forecast  model  (Hogan 
and  Rosmond  1991;  Hogan  et  al.  1991). 
Operational  NOGAPS  forecasts  currently 
consist  of  high-resolution  T239L30  six-day 
forecasts  every  6  hours  and  once-a-day 
extended  ten-day  guidance  using  the  FNMOC 
ensemble  (T119L24).  These  forecasts  are 
distributed  to  numerous  defense  and  civilian 
users,  and  are  used  as  input  for  numerous  DoD 
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environmental  and  application  systems. 
Examples  include:  the  Navy's  Coupled  Ocean- 
Atmosphere  Mesoscale  Prediction  System 
(COAMPS™);  FNMOC's  ocean  wave,  sea 
ice,  ocean  thermodynamics,  and  tropical 
cyclone  models;  aircraft  and  ship-routing 
programs. 

NOGAPS  currently  generates  forecast 
fields  from  the  ground  to  -30  km.  In  October, 
2000,  NRL  initiated  an  applied  research 
project  among  its  Space  Science,  Marine 
Meteorology  and  Remote  Sensing  Divisions  to 
extend  the  vertical  range  of  NOGAPS 
forecasts,  initially  to  -60  km,  then 

progressively  to  -80  km  and  then  to  -100  km. 
This  project  is  ongoing  and  has  led  to  a  new 

prototype  NOGAPS  global  model  with 

Advanced  Level  Physics  and  High  Altitude 
(NOGAPS-ALPHA). 

Developing  NOGAPS-ALPHA  involved 
significant  changes  in  model  configuration 
and  implementation  of  new  algorithms  and 
physics/chemistry  packages  for  the 

stratosphere  and  lower  mesosphere.  This 
paper  summarizes  this  work,  and  also  briefly 
illustrates  some  of  the  new  capabilities  that 
result  (e.g.,  prognostic  ozone).  Other  papers, 
both  in  this  volume  (Allen  et  al.  2004)  and 
elsewhere  (McCormack  et  al.  2004),  will 
focus  more  specifically  on  the  performance  of 
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NOGAPS-ALPHA  as  a  global  forecasting  and 
modeling  tool  for  the  stratosphere. 

2.  THE  NOGAPS  GLOBAL  SPECTRAL 
FORECAST  MODEL  (GSFM) 

Hogan  and  Rosmond  (1991)  and  Hogan  et 
al.  (1991)  provide  thorough  descriptions  of  the 
NOGAPS  forecast  model.  Briefly,  NOGAPS 
is  a  global  spectral  model  (GSM)  in  the 
horizontal  and  utilizes  a  generalized  vertical 
coordinate  within  an  energy  conserving 
vertical  finite  difference  formulation 
(Kasahara  1974;  Simmons  and  Burridge  1981; 
Hogan  and  Rosmond  1991).  The  model’s 
dynamical  variables  are  relative  vorticity, 
divergence,  virtual  potential  temperature, 
specific  humidity,  and  terrain  (surface) 
pressure.  The  model  is  central  in  time  with  a 
semi-implicit  treatment  of  gravity  wave 
propagation  and  Robert  (Asselin)  time 
filtering  (Simmons  et  al.  1978).  The  current 
operational  model’s  physics  packages  include 
a  bulk-Richardson  number  dependent  vertical 
mixing  scheme  (Louis  et  al.  1982),  a  time- 
implicit  Louis  surface  flux  parameterization 
(Louis  1979),  orographic  gravity  wave  drag 
(Webster  et  al.  2003),  shallow  cumulus 
mixing  of  moisture,  temperature,  and  winds 
(Tiedtke  1984),  the  Emanuel  cumulus 
parameterization  (Emanuel  and  Zivkovic- 
Rothman  1999;  Peng  et  al.  2004),  convective, 
stratiform  and  boundary  layer  cloud 
parameterizations  (Slingo  1987;  Teixeira  and 
Hogan  2002),  and  a  shortwave  and  longwave 
radiation  scheme  (Harshvardhan  et  al.  1987). 

3.  NOGAPS-ALPHA  Developments 

Our  current  prototype  NOGAPS-ALPHA 
model  extends  to  ~85  km  in  altitude,  resolving 
the  entire  stratosphere  and  parts  of  the 
mesosphere.  To  accurately  model  these  new 
upper-level  regions  of  the  atmosphere  required 
developing,  implementing  and  testing  of  a 
whole  new  family  of  algorithms  and  physics 


packages  for  the  stratosphere  and  lower 
mesosphere. 

Throughout  this  development  process, 
particular  attention  has  been  paid  to 
maintaining  close  connections  between  the 
NOGAPS-ALPHA  GSM  and  the  continually 
evolving  operational  NOGAPS.  Whenever 
new  tropospheric  physics  packages  have  been 
added  to  NOGAPS  (e.g.,  Hogan  et  al.  2003; 
Peng  et  al.  2004),  these  new  packages  have 
also  been  quickly  transitioned  into  NOGAPS- 
ALPHA.  Additionally,  any  new  physics 
package  added  to  NOGAPS-ALPHA  that 
purports  to  improve  on  an  existing  operational 
scheme  never  results  in  removal  of  the  current 
operational  scheme  from  the  NOGAPS- 
ALPHA  code.  Instead,  NOGAPS-ALPHA 
makes  extensive  use  of  “switches”  which 
enable  all  the  current  operational  and  new 
ALPHA  physics  packages  and  algorithms  to 
be  activated  or  deactivated  in  various  ways  for 
any  given  run.  This  ensures  that  NOGAPS- 
ALPHA  remains  a  high-altitude  prototype  of 
the  current  state-of-the-art  NOGAPS  GSFM 
running  operationally  at  FNMOC. 

We  now  summarize  some  of  the  more 
important  new  features  in  NOGAPS-ALPHA. 

3a.  Hybrid  cr-p  Vertical  Coordinate 

The  NOGAPS  GSM  utilizes  a  vertical 
coordinate  of  the  general  Simmons  and 
Strafing  (1983)  analytical  form 

PClk+iii’ Ps )  =  d(%+1/2)  +  b(rj/c+l/2)[ps  —  Ptop\'  (1*0 

Since  the  normalized  model  coordinate  r/k+i/2 
need  not  be  defined  explicitly  in  (la)  for  use 
in  the  GSM,  we  reexpress  (la)  using  the 
implicit  subscript  notation 

Pk+l/2  =  -4t+l/2  +  Bk+l/2[Ps  PtOP ]’  (Ik) 

where  k=l,..,L  is  the  integer  full  model  layer 
number  (from  top  to  surface),  pk+rn  are  the 
L+l  half  level  interface  pressures,  ps  is  the 
surface  (terrain)  pressure  and  prop  is  the 
(constant)  model  top  pressure.  Through 
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judicious  choices  for  the  coefficients  Ak+m 
and  Bk+1/2  at  each  hybrid  level,  (1)  permits 
hybrid  vertical  model  levels  that  transition 
from  terrain- following  (cr- like)  levels  near  the 
ground  to  pure  pressure  levels  in  the 
stratosphere,  subject  to  the  boundary 
conditions  Pi/2=Ptop  (A  i/2=Ptop,  Bj/2=0, 
Pi/2=0)  and  pL+i/2=Ps  ( Al+i/2=Ptop ,  Bi+m=\, 

f]L+ 1/2=1)  ■ 

NOGAPS  has  always  used  a  coordinate 
defined  in  (lb)  by  setting  Ak+m=Prop  at  all 
levels  k,  so  that5t+1/2=  [p  - pmp]/[ps  - pTOP] , 
(Hogan  and  Rosmond  1991),  yielding  a 
coordinate  almost  identical  to  the  standard  a 
coordinate  p/ps  (since  ps  » pTOP ).  Those 

levels  for  the  current  operational  L30  model 
are  plotted  in  Figure  la. 

In  NOGAPS-ALPHA  we  now  utilize  the 
full  functionality  of  (lb)  by  choosing  Ak+m 
and  Bk+i/2  coefficients  that  give  hybrid  cr-p 
levels  that  transition  smoothly  from  terrain 
following  to  isobaric  in  the  stratosphere:  Fels 
et  al.  (1980),  Simmons  and  Burridge  (1981) 
and  Trenberth  and  Stepaniak  (2002)  provide 


(a)  NOGAPS  OPS  (T239L30)  (b)  NOGAPS-ALPHA  (T239L54) 


34.5°N  for:  (a)  operational  NOGAPS  T239L30 
model,  with  proP=  1  hPa  (red)  and  cr  coordinate;  (b) 
NOGAPS-ALPHA  T239L54  model  with 
Ptop= 0.005  hPa  using  our  new  hybrid  <7-p 
formulation  based  on  (lb)  which  transitions 
smoothly  from  terrain  following  to  a  first  purely 
isobaric  level  at  pINj= 72.6  hPa  (k=  34),  shown  in 
yellow.  The  top  2  full  model  layers  in  NOGAPS 
experience  enhanced  damping  and  diffusion  to  act 
as  a  sponge  layer,  and  the  top  3  model  half  levels 
bounding  these  top  2  layers  are  marked  in  orange. 


motivation  for  using  isobaric  model  levels  in 
the  stratosphere  and  mesosphere. 

Initially  we  implemented  and  tested  a 
Sangster-Arakawa-Lamb  (SAL)  coordinate 
(Sangster  1960;  Arakawa  and  Lamb  1977), 
which  gives  cr-like  levels  up  to  some  interface 
pressure  pINT  (typically  -70-100  hPa)  and  then 
isobaric  levels  from  pINT  to  prop-  The  SAL 
hybrid  coordinate  is  used  in  a  number  of 
global  models:  e.g.,  GFDL  SKYHI  (Fels  et  al. 
1980),  ARPEGE  (Cariolle  and  Deque  1986), 
and  CAM2  (Collins  et  al.  2003). 

However,  the  SAL  coordinate  has  the 
disadvantage  of  producing  sharp 
discontinuities  in  pressure  height  thicknesses 
at  the  pressure  interface  pint  over  high  terrain, 
an  effect  that  can  yield  large  numerical  errors 
in  computed  vertical  gradient  terms  that  are 
used  in  the  governing  equations  (e.g., 
Simmons  and  Burridge  1981). 

Consequently,  for  NOGAPS-ALPHA  we 
have  developed  a  more  elaborate  hybrid 
coordinate  formulation  that  generates  Ak+m 
and  Bk+1/2  coefficients  that  transition  smoothly 
to  isobaric  levels  without  any  abrupt 
discontinuities  in  layer  thicknesses  over 
mountains  (achieved  by  ensuring  that 
{dB/d //) — >0  smoothly  as  p~>Pint)-  The  exact 
equations  and  procedures  used,  though 
straightforward,  are  lengthy  and  will  be 
reported  in  full  elsewhere.  The  hybrid  model 
levels  it  produces  for  our  standard  L54 
NOGAPS-ALPHA  configuration  with 
Ptop= 0.005  hPa  are  plotted  in  Figure  lb, 
illustrating  the  smooth  transition  to  constant 
layer  thicknesses  in  the  stratosphere  regardless 
of  topographic  elevation. 

Operational  NOGAPS  uses  model  layer 
thicknesses  that  dilate  rapidly  with  increasing 
model  height  through  the  stratosphere  (Figure 
la).  By  contrast,  our  working  NOGAPS- 
ALPHA  formulation  utilizes  model  layers 
with  constant  pressure  height  thicknesses  of 
-2  km  throughout  the  stratosphere  and  lower 
mesosphere.  From  Figure  1,  we  see  that 
undamped  model  layers  in  NOGAPS-ALPHA 
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New  Parameterizations  of  Solar  Heating  Rate  of  the  Atmosphere 
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Figure  2:  Profiles  of  diumally-averaged 
shortwave  heating  (SW)  at  45°N  on  21  June 
using:  (a)  a  detailed  middle  atmosphere  SW 
radiation  calculation  from  NRL’s  two- 
dimensional  CHEM2D  model  (blue  curve), 
which  we  use  as  validation  for  two  NOGAPS- 
ALPHA  schemes;  (b)  current  NOGAPS  scheme 
due  to  Harshvardhan  et  al.  (1987)  (black  curve); 
(c)  CLIRAD  SW  scheme  due  to  Chou  and  Suarez 
(2002)  (red  curve).  Note  the  differences  in  the 
upper  stratosphere. 

extend  the  vertical  forecast  range  of  the 
current  operational  model  by  a  factor  of  ~3. 

3b.  CLIRAD  Radiation  Scheme 

To  provide  more  realistic  radiative  driving 
of  the  atmosphere,  we  have  implemented  the 
NASA  “CLIRAD”  longwave  (LW)  and 
shortwave  (SW)  schemes  of  Chou  et  al. 
(2001)  and  Chou  and  Suarez  (2002), 
respectively.  CLIRAD  provides  improved 
radiation  calculations  at  all  model  levels,  but 
significantly  improves  the  middle  atmosphere. 
For  example,  CLIRAD  SW  heating  rates 
include  contributions  from  O2  and  near- 
infrared  CO2  bands  that  are  not  contained  in 
the  current  operational  radiation  scheme 
(Harshvardhan  et  al.,  1987),  which  lead  to 
underestimates  in  peak  SW  heating  in  the 
middle  atmosphere  of  as  much  as  1-2  K  day'1, 
as  shown  in  Figure  2. 

To  specify  the  global  SW  heating  and  LW 
cooling  rate  contributions  from  ozone  (O3) 
and  water  vapor  (H20),  NOGAPS  utilizes 


seasonally-varying  zonal-mean  pressure- 
latitude  lookup  tables  of  their  climatological 
mixing  ratios,  continually  interpolating  these 
values  to  the  current  forecast  time  and  onto 
the  NOGAPS  model  levels  and  Gaussian  grid. 
For  NOGAPS-ALPHA,  the  existing  NOGAPS 
O3  climatology  of  McPeters  et  al.  (1984)  was 
updated  and  extended  vertically  using  the 
monthly  mean  O3  climatology  of  Fortuin  and 
Kelder  (1998)  (1000-0.3  hPa),  while  upper- 
level  water  vapor  was  prescribed  using 
measurements  (100-0.3  hPa)  from  the 
Halogen  Occultation  Experiment  (HALOE) 
(Randel  et  al.  2001).  At  heights  from  0.3  hPa 
up  to  0.005  hPa,  both  O3  and  H2O  mixing 
ratios  are  based  on  long-term  climate  output 
from  NRL’s  zonally-averaged  CHEM2D 
dynamics-radiation-chemistry  global  model, 
which  spans  the  -0-108  km  height  range  and 
includes  very  detailed  middle  atmospheric 
chemistry  and  radiation  packages  (Summers  et 
al.  1997;  McCormack  and  Siskind  2002). 

Our  goal  is  to  ultimately  replace  these  2D 
climatological  mixing  ratios  in  the  radiation 
calculations  with  three-dimensional  (3D) 
NOGAPS-ALPHA  prognostic  chemical 
constituent  fields,  which  are  introduced  and 
discussed  in  section  4.  In  NOGAPS-ALPHA 
we  can  already  choose  to  use  our 
developmental  3D  prognostic  ozone  fields 
instead  of  2D  climatology  within  the  CLIRAD 
radiation  scheme  to  specify  the  O3  LW  and 
SW  contributions. 

3c.  Upper  Level  Meteorological  Initialization 

In  2003,  FNMOC  replaced  its  Multivariate 
Optimal  Interpolation  (MVOI)  data 
assimilation  system  (Goerrs  and  Phoebus 
1992)  with  the  new  NRL  Atmospheric 
Variational  Data  Assimilation  System 
(NAVDAS),  which  is  based  on  three- 
dimensional  variational  (3DVAR)  analysis 
(Daley  and  Barker  2001).  MVOI  fields 
extended  only  to  10  hPa,  though  an 
experimental  “STRATOI”  product  was  also 
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issued,  based  mostly  on  TIROS-N  Operational 
Vertical  Sounder  (TOVS)  radiances,  which 
yielded  additional  initialization  wind, 
temperature  and  geopotential  fields  up  to  0.4 
hPa.  Our  test  runs  with  NOGAPS-ALPHA  to 
date  have  used  MVOI/STRATOI  fields  almost 
exclusively.  Newer  NAVDAS  fields  currently 
extend  to  ~4  hPa,  but  should  later  extend  to 
~0.1  hPa  when  direct  radiance  assimilation 
comes  online.  Thus,  to  initialize  NOGAPS- 
ALPHA  we  must  for  now  supplement  current 
operational  meteorological  initialization  fields 
with  some  working  representation  of  the 
atmospheric  state  between  ptop=0.005  hPa 
and  the  top  operational  analysis  level  pan- 

Our  current  approach  is  to  extrapolate  the 
topmost  analysis  temperatures  and  zonal  and 
meridional  wind  fields  upwards  and  then 
progressively  relax  the  extrapolated  profiles 
with  increasing  altitude  towards  monthly 
varying  zonal-mean  climatological  values.  For 
temperature,  the  equations  used  are 

T(A(p,P)  =  [1  -  H<  p)\  Tih  ( 1,  (p,  p)  +  vt(  pWcum  ( <P,  P) 

(2a) 

V.  (2b> 

1  Pa\  In  pTOP  J 

where  A  is  longitude  and  rp  is  latitude. 
Tclim(<P,p)  are  zonal-mean  climatological 
temperatures.  TLE(A,(p,p)  are  “lapse-rate 
extrapolated”  temperatures  in  which  the 
analysis  temperatures  at  the  top  analysis  level, 
T(X,(p,pan),  are  extrapolated  upwards  using  the 
climatological  lapse  rates  dTcLm/dp.  Figure  3a 
shows  how  the  final  profile  T(A,rp,p)  at 
altitudes  above  pan  is  formed  as  the  weighted 
sum  of  these  two  profiles,  as  specified  in  (2). 

The  normalized  weighting  coefficient  w(p) 
in  (2b)  increases  monotonically  with 
increasing  log-pressure  altitude,  such  that 
w(Pan)= 0  and  w(prop)= 1,  so  that  the  fit  (2a) 
contains  pure  analysis  at  pan,  mostly 
extrapolated  analysis  near  pan  (to  retain  close 
connections  with  analyzed  fields),  but  relaxes 


to  pure  climatology  at  prop-  The  rate  at  which 
T(A,rp,p)  reverts  to  climatology  is  controlled 
by  the  coefficient  s  in  (2b).  For  small  (large) 
s,  (2a)  transitions  more  rapidly  (gradually) 
with  altitude  to  climatology.  Currently  we  use 
e=\  in  NOGAPS-ALPHA.  The  climatological 
fields  we  use  are  monthly  zonal-mean  zonal 
winds,  temperatures  and  geopotential  heights 
from  the  1986  COSPAR  International 
Reference  Atmosphere  (CIRA)  (Fleming  et  al. 
1990),  interpolated  to  the  initialization  time 
and  model  pressures  and  latitudes.  The  final 
temperature  profile  (2a)  in  turn  yields  a 
geopotential  profile  though  the  hydrostatic 
relation 

</>(X,(p,p)  =  <j)AN+R\PT0P  T  (A, (p,p)d\np  (3) 

—  Pan 

where  (f>AN=(f{X,(p,pAN)  and  R  is  the  gas  constant. 

Winds  are  dealt  with  slightly  differently  to 
temperatures.  As  shown  in  Figure  3b,  we 
perform  a  direct  upward  extrapolation  of 

20  December  1999  00:00Z:  40.5' W,  69.3  N 


(a)  Initialization  Temperature  T(p)  (b)  Initialization  Zonal  Wind  u(p) 


Figure  3:  Black  curve  shows  NOGAPS-ALPHA 
initialization  of  (a)  temperature  and  (b)  zonal 
winds  for  20  December  1999  0:00Z  at  a  sample 
location  of  40.5°W,  69.5°N.  Other  curves  show: 
NOGAPS  MVOI/STRATOI  operational  analysis 
(aqua  curve);  profiles  extrapolated  upwards  from 
top  analysis  level  using  CIRA  lapse  rates  (dark 
green)  and  constant  wind  (light  green).  Purple 
curves  show  CIRA  climatology  for  December  at 
69.5°N.  Final  profile  is  a  progressively  weighted 
fit  between  the  green  and  purple  curves,  yielding 
pure  analysis  (green)  at  0.4  hPa  and  pure  CIRA 
climatology  (purple)  at  prop= 0.005  hPa. 
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winds  at  the  top  analysis  levels  and 
progressively  blend  these  with  climatological 
winds  using  (2b).  Given  the  unreliability  of 
mesospheric  CIRA  zonal  winds  during  certain 
months,  particularly  near  the  equator  (Randel 
et  al.  2004),  we  have  added  an  additional 
option  of  relaxing  zonal  winds  to  either  CIRA 
or  the  more  recent  climatological  zonal  wind 
fields  from  the  UARS  Reference  Atmosphere 
Project  (URAP)  (Swinbank  and  Ortland 
2003).  Since  we  have  no  reliable  meridional 
wind  climatology  at  upper  altitudes,  we 
currently  relax  meridional  winds  to  a 
reference  climatology  field  set  to  zero 
everywhere,  yielding  vanishing  meridional 
winds  initially  at  prop- 

In  addition  to  our  standard  operational 
NAVDAS  fields,  we  have  developed 
strategies  for  ingesting  into  NOGAPS- 
ALPHA  initialization  fields  from  other  NWP 
centers,  such  as  the  European  Centre  for 
Medium-Range  Weather  Forecasts  (ECMWF) 
and  the  National  Centers  for  Environmental 
Prediction  (NCEP):  see  Allen  et  al.  (2004)  for 
more  details. 

3d.  Gravity  Wave  Drag 

The  importance  of  parameterized  subgrid- 
scale  orographic  gravity  wave  drag  (GWD) 
for  improved  forecast  skill  in  the  troposphere 
is  well  known  (e.g.  Kim  et  al.  2003;  Webster 
et  al.  2003).  On  extending  the  vertical  range  of 
NWP  models  into  the  stratosphere  and 
mesosphere,  GWD  parameterization  becomes 
an  even  more  important  and  general  issue, 
since  the  circulations  of  the  entire  equatorial 
stratosphere  and  global  mesosphere  are  driven 
by  large  synoptic-scale  body  forces 
maintained  by  quasi-continuous  gravity  wave 
breaking  around  the  globe  (Fritts  and 
Alexander  2003).  Even  high-resolution  global 
NWP  models  like  NOGAPS  cannot 
adequately  resolve  this  GWD  (Kim  et  al. 
2003). 


For  NOGAPS-ALPHA  we  have  developed 
and  implemented  a  number  of  new  GWD 
schemes.  We  have  implemented  three  separate 
orographic  GWD  schemes,  any  one  of  which 
can  be  activated  for  a  given  run.  Two  are 
operational  schemes:  the  Palmer  et  al.  scheme 
(Palmer  et  al.,  1986),  which  was  used 
operationally  in  NOGAPS  for  many  years,  and 
the  newer  scheme  of  Webster  et  al.  (2003), 
which  replaced  it  in  2003  (Hogan  et  al.,  2003). 
The  third  scheme  is  a  more  complex 
developmental  scheme  of  Kim  and  Doyle 
(2004)  that  extends  earlier  work  of  Kim  and 
Arakawa  (1995),  a  version  of  which  is 
implemented  currently  in  the  NCEP  GSM 
(Alpert  et  al.  1996).  Kim  and  Hogan  (2004) 
report  on  momentum  budget  studies  with 
NOGAPS  using  this  scheme.  A 

parameterization  of  GWD  due  to  convectively 
generated  waves  has  also  been  developed  and 
installed  (Chun  and  Baik  2002). 

In  the  stratosphere  and  mesosphere, 
gravity  wave  amplitudes  become  large  and 
there  is  greater  interaction  among  a  larger 
number  of  waves,  leading  to  less  obvious 
signatures  of  individual  wave  packets  from 
identifiable  tropospheric  sources  like 

mountains  and  deep  convection  (Fritts  and 
Alexander,  2003).  Thus,  we  have  also 

implemented  and  begun  testing  several 

spectral  GWD  schemes  that  have  become 
popular  recently  for  parameterizing  GWD 
from  a  broad  spectrum  of  waves  in  the  middle 
stratosphere  and  above.  We  have  implemented 
the  Alexander  and  Dunkerton  (1999)  scheme 
into  NOGAPS-ALPHA,  generalizing  it  to 
multiple  wave  propagation  azimuths.  We  have 
tested  it  in  a  series  of  12-24  month  runs  in 
efforts  to  tune  it  to  generate  a  realistic 
semiannual  zonal  wind  oscillation  (SAO)  in 
the  tropical  stratosphere  and  lower 
mesosphere. 

We  have  also  investigated  coupling  among 
various  schemes  that  might  provide  a  more 
integrated  GWD  product  for  NOGAPS- 
ALPHA.  We  have  coded  the  so-called 
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“Doppler-spread”  spectral  GWD  scheme  of 
Hines  (in  collaboration  with  B.  N.  Lawrence 
from  Rutherford-Appleton  Laboratory)  and 
can  couple  it  with  our  orographic  GWD 
schemes,  in  the  general  manner  outlined  by 
Hines  (1997).  We  have  also  investigated  ways 
in  which  our  orographic  and  convective  GWD 
schemes  lower  down  might  be  coupled  to  an 
Alexander-Dunkerton  spectral  scheme  higher 
up.  However,  given  the  complexity  of  these 
developmental  coupled  schemes,  we  are 
currently  testing  them  offline  in  single-column 
tests  using  NOGAPS-ALPHA  model  grids 
and  fields. 

Figure  4  shows  an  example  of  single 
column  NOGAPS-ALPHA  tests  of  some  of 
these  individual  (uncoupled)  GWD  schemes. 
This  case  uses  model  fields  at  +114  hours 
from  a  NOGAPS-ALPHA  T79L54  forecast 
run  initialized  from  MVOI/STRATOI  fields 
on  17  January  2003  at  0:00Z.  Figures  4a  and 
4b  show  zonal  wind  and  temperature, 
respectively,  over  southern  Scandinavia.  At 
this  time  (21  Feburary,  18:00Z)  a  dynamical 
stratospheric  warming  event  was  well 
underway  and  the  vortex  temporarily  split 
(e.g.,  deZaffa  and  Muscari  2004;  McCormack 
et  al.  2004).  The  warming  is  evident  from  the 
enhanced  NOGAPS-ALPHA  stratopause 
temperature  compared  to  CIRA  climatology. 

The  middle  row  of  plots  shows  mean-flow 
accelerations  (MFAs)  due  to  parameterized 
orographic  drag  from  flow  across  the  southern 
Kjonas  Mountains,  as  specified  by  the  Palmer 
et  al.  GWD  scheme  (Figure  4c)  and  the 
Webster  et  al.  scheme  (Figure  4d).  The 
differences  are  broadly  consistent  with  some 
large  differences  in  physics  between  these  two 
schemes:  the  newer  Webster  et  al.  scheme 
produces  much  larger  flow  blocking  drag  near 
the  surface  and  transfers  a  much  smaller 
fraction  of  the  total  surface  pressure  drag 
across  the  mountains  into  propagating 
mountain  waves,  compared  to  the  Palmer  et  al. 
scheme  (Webster  et  al.,  2003).  This  explains 
the  significant  surface  MFA  and  much  smaller 


NOGAPS-ALPHA  T79L54  Forecast:  Initialized  17  December  2003  0:00Z 
+114  Hour  Forecast:  21  December  2003  18:00Z:  Profile  at  7.5°E,  60.5°N 
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Figure  4:  (a)  Zonal  wind  and  (b)  temperature 
over  southern  Scandinavia  (7.5°E,  60.5°N)  from  a 
NOGAPS-ALPHA  T79L54  +114  hour  forecast 
initialized  on  17  January  2003  at  0:00Z.  Middle 
row  shows  zonal  mean-flow  accelerations  from 
orographic  GWD  from  (c)  Palmer  et  al.  and  (d) 
Webster  et  al.  schemes.  Bottom  row  shows  mean- 
flow  accelerations  from  spectral  GWD  schemes 
of  (e)  Alexander-Dunkerton  (AD)  and  (f)  Hines. 
The  AD  profiles  show  sensitivity  of  profile  to 
number  of  phase  speeds  nc  in  the  spectmm.  Red 
curve  shows  pTop  and  yellow  curve  shows  pINT- 


stratospheric  MFAs  in  Figure  4d  compared  to 
the  Palmer  et  al.  profile  in  Figure  4c. 

Figure  4e  shows  a  series  of  MFA  profiles 
computed  using  the  spectral  scheme  of 
Alexander  and  Dunkerton  (1999),  illustrating 
one  of  numerous  tuning  experiments  needed 
when  implementing  such  schemes.  These 
profiles  investigate  the  sensitivity  of  the  MFA 
profile  to  one  particular  parameter:  the 
number  of  individual  phase  speeds  nc  in  the 
wave  spectrum.  We  see  the  profiles  generally 
overlay  for  nc^.  120.  This  wave  spectrum  was 
launched  at  15  km,  with  a  phase  speed  width 
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cw=60  m  s'1  centered  about  the  local  mean 
wind  at  15  km,  peak  momentum  flux  per  unit 
mass  Bm=0A  m  s'  ,  a  horizontal  wavelength  of 
100  km,  and  total  flux  per  unit  mass 
F50=0.0004m  s'  .  This  spectrum  is  launched  at 
3  different  propagation  azimuth  pairs:  0/180°, 
45/225°,  and  -45/135°.  All  these  parameters 
must  be  tuned  carefully  too  for  acceptable 
middle  atmospheric  circulations  (e.g., 
Alexander  and  Rosenlof,  2003). 

Figure  4f  illustrates  a  similar  run  with  the 
Hines  (1997)  spectral  scheme.  The  MFA 
profile  is  generally  similar  though  less  smooth 
vertically,  consistent  with  findings  in  some 
previous  studies  which  show  that  the  Hines 
scheme  deposits  wave  fluxes  more  abruptly  at 
high  altitudes  rather  than  smoothly  over  a 
more  extended  altitude  range  (Charron  et  al. 
2002). 

As  with  lower  atmospheric  schemes,  these 
upper  atmospheric  GWD  schemes  need  to  be 
carefully  tuned  to  yield  acceptable  circulations 
(e.g.,  Figure  4e),  which,  given  their 
complexity,  is  a  time-consuming  process  and 
necessitates  retuning  whenever  horizontal  or 
vertical  model  resolution  is  changed.  Rayleigh 
friction  (RF)  can  also  be  used:  while  RF  is  a 
much  cruder  proxy  for  mesospheric  GWD  and 
can  have  deleterious  side  effects,  such  as 
suppression  of  local  dynamical  variability  and 
the  potential  for  unrealistic  secondary 
circulations  (Shepherd  et  al.  1996),  it  has  the 
advantage  of  being  more  easily  tunable  and  so 
is  used  in  the  upper  boundaries  of  many  global 
NWP  models  (e.g.,  Butchart  and  Austin 
1998).  Thus  we  have  also  implemented  a 
general  RF  scheme  for  simpler  middle- 
atmospheric  GWD  specification,  which 
currently  allows  the  user  to  choose  from  any 
one  of  nine  prescribed  RF  profiles.  As  shown 
in  Figure  5,  these  profiles  range  from  weak  to 
strong  and  (vertically)  shallow  to  deep  drag, 
and  are  all  based  on  profiles  used  in  other 
global  middle-atmosphere  model  studies 
(Boville  1986;  Shepherd  et  al.  1996;  Butchart 


Rayleigh  Friction  Profiles 


1.  Boville  (J.  Atmos.  Sci.,  1986) - 

2.  Kim-modified  Boville  (1986) - 

3.  Lin  and  Williamson  (2000) - 

4.  Shepherd  et  al.  (J.  Geophys.  Res.,  1996) - 

5.  Standard  Butchart  and  Austin  (J.  Atmos.  Sci.,  1998) - 

7.  Limpasuvan  et  al.  (J.  Atmos.  Sci.,  2000) - 

8.  Normal  Pawson  et  al.  (Q.  J.  R.  Meteorol.  Soc.,  1998) - 

9.  Deep  Pawson  et  al.  (Q.  J.  R.  Meteorol.  Soc.,  1998) - 
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Figure  5:  Various  Rayleigh  friction  profiles 
(expressed  as  damping  rates)  available  for  use  in 
NOGAPS-ALPHA,  based  on  profiles  used  in 
some  other  global  models  and  spanning  various 
ranges  of  required  damping.  Each  profile  was 
chosen  for  a  specific  model  with  a  differing  prop, 
and  those  tops  for  each  model  are  marked  by  a 
change  from  solid  to  dotted  curves  in  each  case. 


and  Austin  1998;  Pawson  et  al.  1998; 
Limpasuvan  et  al.  2000). 


4.  FULL  PROGNOSTIC  CAPABILITY  FOR 
TRACE  CHEMICAL  CONSTITUENTS 


Another  major  addition  to  NOGAPS- 
ALPHA  has  been  full  3D  prognostic 
capabilities  for  key  trace  chemical 
constituents,  such  as  ozone.  Just  as  for 
meteorological  variables,  full  prognostic 
chemistry  capability  requires  many  features, 
including  3D  global  initialization,  accurate 
global  advection  (adiabatic  tendencies)  and, 
where  necessary,  photochemical  and/or 
microphysical  updating  (diabatic  tendencies). 
We  devote  a  separate  section  here  to 
describing  this  new  capability. 

4a.  Initialization 


The  number  of  specific  chemicals  to  be 
used  in  any  given  NOGAPS-ALPHA  run  is 
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specified  by  various  logical  switches  at  run 
time.  NOGAPS-ALPHA  currently  supports 
initialization  (and  subsequent  forecasting)  of 
the  following  minor  constituents:  water  vapor 
(H2O),  ozone  (O3),  nitrous  oxide  (N2O),  and 
methane  (CH4).  In  addition,  each  chemical  can 
be  forwarded  in  either  a  “passive”  mode  (no 
chemistry/microphysics,  purely  passive 
transport),  or  “active”  mode 
(chemistry/microphysics  activated),  allowing 
for  a  total  of  8  different  types  of  prognostic 
chemical  fields.  The  passive  and  active 
chemical  pairs  prove  particularly  valuable 
when  assessing  new  parameterized  chemistry 
schemes  in  NOGAPS-ALPHA,  since 
subtracting  the  two  fields  allows  us  to  factor 
out  transport  and  thus  isolate  the  accumulated 
chemistry  effects  within  air  parcels  at  specific 
locations. 

However,  active  water  vapor  has  always 
been  an  intrinsic  model  variable  in  NOGAPS, 
since  it  enters  the  GSM’s  governing  equations 
of  motion  through  the  virtual  potential 
temperature  (Hogan  and  Rosmond  1991): 
thus,  the  model  treats  “active  H2O”  as  a 
meteorological  field,  reducing  the  total 
number  of  possible  chemical  fields  to  7, 
including  passive  H20. 

Water  vapor  is  initialized  as  in  NOGAPS, 
using  3D  analyzed  dew  point  depression 
fields,  which  are  subsequently  converted  to 
specific  humidity  fields.  These  fields  do  not 
currently  extend  much  above  10  hPa  and  so 
the  10  hPa  analyzed  humidity  fields  are 
simply  extrapolated  to  upper  altitudes  in 
NOGAPS-ALPHA,  which,  given  the  relative 
dryness  of  the  stratosphere  and  mesosphere,  is 
an  acceptable  short-term  solution. 

Ozone  can  be  initialized  in  one  of  two 
ways  in  NOGAPS-ALPHA.  While  NAVDAS 
does  not  at  present  assimilate  O3 
operationally,  there  are  plans  to  do  and  so  we 
have  developed  code  in  NOGAPS-ALPHA  to 
ingest  and  initialize  such  3D  initialization 
fields  whenever  they  become  available.  In  the 
meantime,  to  test  and  utilize  this  capability  we 


have  generated  proxy  three-dimensional 
analyzed  O3  mixing  ratio  fields  for  NOGAPS- 
ALPHA  runs  based  on  the  O3  analyses  issued 
by  both  the  NASA  Global  Modeling  and 
Analysis  Office  (GMAO)  (Stajner  et  al.  2001) 
and  the  ECMWF  (Dethof  2003).  When 
activated,  NOGAPS-ALPHA  searches  for  3D 
ozone  pressure-level  initialization  fields  for 
the  date  in  question,  reads  them  in  and  then 
interpolates  them  onto  model  levels  and  the 
Gaussian  grid.  Since  these  analyzed  ozone 
fields  currently  abate  at  pan  ~  0.2-1  hPa,  for 
now  we  extrapolate  O3  mixing  ratios  at  this 
top  analysis  level  upwards  to  prop  in  a  manner 
similar  to  zonal  wind,  as  shown  in  Figure  3b 
(green  curve).  In  future  we  will  develop  a 
better  extrapolation  procedure  that  relaxes  to 
the  2D  climatology,  in  a  manner  somewhat 
similar  to  that  for  temperature  fields,  as 
illustrated  in  Figure  3  a. 

When  3D  O3  analysis  fields  are  not 
available,  we  can  instead  initialize  our  3D 
chemical  field  using  the  zonally  averaged  two- 
dimensional  (2D)  NOGAPS-ALPHA  03 
climatological  fields  described  in  section  3b. 

For  N2O  and  CH4,  operational  3D  analyses 
do  not  currently  exist,  and  so  our  current 
strategy  is  to  initialize  these  fields  using  2D 
climatologies.  Our  2D  N2O  and  CH4 
climatologies  for  NOGAPS-ALPHA  were 
derived  from  multi-year  simulations  using 
CHEM2D  (McCormack  and  Siskind  2002), 
and  extend  from  1000  to  0.001  hPa. 

4b.  Chemical  Advection 

Once  initialized,  these  chemical  fields  are 
transported  globally  within  the  GSM  using  the 
same  spectral  advection  algorithms  used  to 
advect  meteorological  variables  such  as 
vorticity  and  specific  humidity.  Over  time, 
spectral  advection  yields  fine  scale  spectral 
noise  in  the  chemical  fields  (Rasch  and 
Williamson  1990;  Fairlie  et  al.  1994),  which 
we  effectively  suppress  through  the  use  of  a 
small  amount  of  V4  horizontal  spectral 
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3  January  2002  0:00Z  Active  N20  (ppbv)  A/nriPC  fv 

Forecast  Hour  +000336  Pressure  8.669  hPa  (Altitude  -33.2  km)  ^  UU/\r  O  -  CX 


Figure  6:  NOGAPS-ALPHA  +336  hour  forecasts 
of  active  N20  at  8.7  hPa  based  a  run  initialized  on 
20  December  2001  at  0:00Z.  Top  plot  shows 
T79L54  run,  bottom  plot  shows  T239L54  run. 
Values  in  color  bar  are  mixing  ratios  in  ppbv, 
contours  are  geopotential  height  in  kilometers 

diffusion  after  each  time  step.  This  can  be 
applied  either  to  the  isobaric  stratospheric 
levels  only,  or  to  all  the  hybrid  model  levels. 

4c.  Photochemistry 

For  direct  simulations  of  stratospheric 
ozone  photochemistry,  a  reasonable 
simulation  would  require  at  least  60  different 
chemical  reactions  involving  ~40  different 
minor  chemical  species.  At  this  early  stage  of 
our  development,  this  is  impractical  and  so  we 
have  pursued  simpler  parameterization 
strategies  that  yield  acceptably  accurate 
photochemical  tendencies  over  forecast  time 
scales. 


We  have  focused  to  date  on  linearized 
ozone  photochemistry  schemes  that 
parameterize  the  major  production/loss  rate 
dependences  globally  through  climatological 
lookup  tables,  which  are  generated  by  a 
middle  atmospheric  climate  model  containing 
sophisticated  radiation  and  chemistry 
packages.  For  NOGAPS-ALPHA,  we  use 
NRL’s  CHEM2D  model  for  this  purpose 
(Summers  et  al.  1997;  McCormack  and 
Siskind  2002). 

Nitrous  Oxide 

We  implemented  prognostic  N20  in 
NOGAPS-ALPHA  primarily  because  its 
photochemistry  is  relatively  simple,  and  thus  it 
is  a  good  candidate  chemical  for  developing 
an  initial  parameterized  photochemistry 
scheme.  However,  prognostic  N2O  is  a  useful 
quantity  to  implement  for  other  reasons  too:  it 
has  important  radiative  effects,  destroys  ozone 
(indirectly,  via  NOx  products  from  N2O 
photodissociation),  and  has  slow 
photochemical  rates  and  sharp  gradients  both 
vertically  and  latitudinally  that  make  it  an 
excellent  diagnostic  pseudo-tracer  of  the 
atmospheric  motion.  Furthermore,  in  the 
stratosphere  it  is  well  measured  by  a  number 
of  satellite  instruments  (e.g.,  Urban  et  al. 
2004)  and  thus  its  prognostic  output  from 
NOGAPS-ALPHA  runs  can  be  validated 
globally. 

N2O  is  produced  at  the  surface  by 
biological  and  anthropogenic  activity  (e.g., 
Huttunen  et  al.  2003;  Scanlon  and  Kiely  2003) 
and  is  destroyed  in  the  stratosphere  primarily 
by  photolysis  (Brasseur  and  Solomon  1984). 
In  parameterizing  this  loss,  we  follow  the 
approach  of  Randel  et  al.  (1994)  and 
parameterize  it  using  the  rate  equation 

d  Xn2o  _  _  Xn2o  ^ 

dt  TNi0((p,p,t) 

where  d/dt  is  the  advective  time  derivative, 
XNi0(X<P,pd)  are  N20  mixing  ratios,  and 
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tno (jp,p,t)  are  the  N2O  loss  rates:  the  latter 

are  derived  from  diumally-averaged  output 
from  long-term  CHEM2D  runs,  then  output  as 
lookup  tables  as  a  function  of  latitude, 
pressure  and  month.  These  loss  rates  are 
generally  slow,  peaking  in  the  equatorial 
stratosphere  at  altitudes  -30-40  km.  The 
tabulated  loss  rate  values  are  continually 
interpolated  within  NOGAPS-ALPHA  to  the 
precise  location  and  forecast  time  and  date  in 
question.  We  approximate  source  production 
by  keeping  the  surface-level  mixing  ratios  set 
to  320  ppbv. 

Figure  6  shows  prognostic  N2O  fields 
after  14  days  in  two  NOGAPS-ALPHA  runs 
at  T79L54  and  T239L54,  with  otherwise 
identical  initial  conditions  and  model 
parameter  settings.  Both  plots  show  a  strong 
wave  breaking  event  in  which  high-N20 
tropical  air,  typically  isolated  from  mid¬ 
latitudes  by  a  subtropical  transport  barrier 
(Neu  et  al.  2003),  is  rapidly  pulled  out  to 
midlatitudes  in  a  thin  streamer  that  then  wraps 
around  into  a  vortex-like  structure  over 
western  Europe.  Close  inspection  of  both  plots 
shows  considerable  similarity,  particularly  in 
the  Southern  Hemisphere,  illustrating  the 
robust  nature  of  these  prognostic  N2O  fields, 
and  in  particular  the  fidelity  of  the  NOGAPS- 
ALPHA  spectral  code  for  chemistry 
advection.  Slight  differences  are  apparent:  the 
streamer  event  in  the  T239L54  run  is 
longitudinally  displaced  from  a  similar 
looking  event  in  the  T79L54  run,  which  is  a 
meteorological  effect  possibly  associated  with 
greater  explicitly  resolved  GWD  in  the 
T239L54  run  that  slightly  modifies  the  mean- 
flow  evolution  between  the  two  runs. 

Methane 

Our  success  in  parameterizing  N2O 
photochemistry  motivated  us  to  develop  a 
similar  chemistry  scheme  for  CH4.  Like  N2O, 
CH4  is  radiatively  important,  affects  odd- 
oxygen  photochemistry,  and  is  measured 


3  January  2002  0:00Z  Active  CH4  (ppbv)  Af/Or/1  PC  rv 

Forecast  Hour  +000336  Pressure  8.669  hPa  (Altitude  -33.2  km)  ^  (jAr  O  -  (X 


0.00000  1500.00 


Figure  7:  NOGAPS-ALPHA  +336  hour  forecasts 
of  active  CH4  at  8.7  hPa  based  a  T79L54  ran 
initialized  on  20  December  2001  at  0:00Z.  Values 
in  color  bar  are  mixing  ratios  in  ppbv,  contours 
are  geopotential  height  in  kilometers. 

globally  from  space.  Additionally,  methane 
oxidation  is  an  important  source  of  water 
vapor  for  the  upper  stratosphere:  a  simple  one¬ 
dimensional  methane  oxidation 

parameterizaton  was  instituted  in  the  ECMWF 
model  to  account  for  this  high-altitude 
moisture  source  for  their  3D  prognostic  water 
vapor  fields.  Thus,  an  ability  to  generate 
prognostic  CH4  in  NOGAPS-ALPHA  could 
then  also  allow  this  stratospheric  moisture 
source  to  be  parameterized  three- 
dimensionally  through  conservation  of  the 
“potential  hydrogen”  term  2  %Ch4+Xh2o 

(Summers  et  al.  1997). 

Similar  to  N2O,  CH4  is  produced  at  the 
surface  and  is  destroyed  in  the  middle 
atmosphere  primarily  by  reactions  with  OH, 
Ot'D),  chlorine  and  Lyman-a  photolysis 
(Brasseur  and  Solomon  1984).  We  have 
parameterized  CH4  chemistry  in  NOGAPS- 
ALPHA  using  an  equation  identical  to  form  to 
(4),  with  the  zCHJ(p,p,t)  diumally  averaged 

loss  rates  again  generated  from  CHEM2D 
runs.  Figure  7  plots  prognostic  CH4  mixing 
ratios  after  14  days  for  a  T79L54  NOGAPS- 
ALPHA  run  with  the  same  meteorological 
initial  conditions  and  very  similar 
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parameterization  settings  to  the  N2O  T79L54 
run  depicted  in  Figure  6,  and  shows  a  similar 
streaming  of  high  CH4  tropical  air  to 
midlatitudes  by  a  wave  breaking  event. 

Ozone 

Ozone  photochemistry  requires  a 
considerably  more  complicated 

parameterization,  since  it  is  continuously 
created  and  destroyed  in  the  middle 
atmosphere.  In  our  initial  implementations,  we 
have  pursued  a  linearized  photochemistry 
approach,  which  begins  by  postulating  a  set  of 
core  model  variables  that  have  the  greatest 
perceived  impact  on  net  ozone  (odd-oxygen) 
production  and  loss  rates,  (P-L).  Extensive 
stratospheric  ozone  research  has  established 
that  production/loss  of  ozone  is  most  sensitive 
to  the  current  local  value  of  the  ozone  mixing 
ratio,  %03,  the  local  temperature  T,  and  the 
local  intensity  of  UV  radiation,  which  in  turn 
depends  on  the  overhead  vertical  column 
amount  of  ozone,  C03,  via  the  shielding  effect. 
We  express  this  functionally  as 

^  =  (P-L)[XvT,c0i]  .  (5) 


3  partial  derivative  terms  in  (6)  can  all  be 
evaluated  with  global  models  containing  full 
radiation  and  chemistry  (e.g.,  CHEM2D)  via  a 
series  of  simulations  that  evaluate  net  rates 
about  this  reference  state  and  then  perturb  that 
state  with  respect  to  the  mixing  ratio, 
temperature,  and  column  amount  to  estimate 
the  partial  derivative  terms.  The  4  linearized 
photochemical  coefficients  in  (6)  can  then  be 
progressively  estimated  and  issued  as  latitude- 
pressure-month  lookup  tables. 

Observational  validation  for  this  linearized 
approach  (6)  to  parameterizing  O3 
photochemistry  exists.  Linearized 

photochemistry  has  been  used  successfully  to 
model  observed  ozone-temperature 
correlations  in  the  stratosphere  (e.g.,  Douglass 
et  al.  1985;  Stolarski  and  Douglass  1985)  and 
various  types  of  stratospheric  planetary  wave 
signals  seen  in  satellite  ozone  data  (e.g., 
Randel  and  Gille  1991;  Randel  1993;  Smith  et 
al.  2002). 

We  have  developed  and  implemented  3 
separate  O3  photochemistry  schemes  for 
NOGAPS-ALPHA  based  on  this  particular 
parameterization  approach,  any  one  of  which 
can  be  selected  by  the  user  for  prognostic 
ozone  runs. 


We  approximate  the  function  on  the  right- 
hand  side  of  (5)  using  a  standard  truncated 
(linearized)  Taylor  series  expansion 
(McLinden  et  al.  2000) 


(P-L)  =  (P-L\  + 


djP-L) 

dZo3 


(Zo3-Zo3)  + 


d( P-L) 


dT 


(t-t)+^p  ^ 


dc, 


o, 


(cQ-coJ-(6) 


The  subscript  “0”  denotes  values  computed  at 
some  reference  state  %0  >  T  ,  and  c0j ,  which 

ideally  should  be  a  chemical  equilibrium  state 
but,  when  implemented  as  a  parameterization 
in  models,  tends  to  work  better  using 
reference  climatologies.  Thus  ( P-L)  and  its 


CHEM2D  scheme 

As  with  our  N2O  and  CH4  schemes,  we 
have  run  CHEM2D  to  generate  O3  linearized 
photochemical  coefficients.  An  important 
point  to  note  is  that,  while  our  focus  is  on 
ozone,  the  linearized  photochemical 
coefficients  themselves  relate  most  directly  to 
the  odd-oxygen  (Ox)  photochemistry  (O3  +  O) 
that  governs  net  O3  production/loss,  rather 
than  the  individual  O3  rates  themselves.  At  an 
equilibrium  Ox  state,  there  is  still  constant 
rapid  chemical  cycling  between  ozone  and 
monatomic  oxygen,  so  that  the  lifetime  of 
individual  O3  molecules  is  short  yet  total  O3 
mixing  ratios  do  not  change  (Brasseur  and 
Solomon  1984).  The  CHEM2D  runs  needed  to 
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generate  these  linearized  O3  photochemical 
coefficients  are  much  more  challenging  and 
time-consuming  than  for  N2O  and  CH4,  and  to 
date  we  have  generated  and  validated  only  the 
first  two  linearized  photochemistry  terms  in 
(6).  We  express  these  here  as 


dXo, 

dt 


/  r>  t\CHEM2D\  „  .1 

(P-L)  0  [<P,P,t\ 


CHEM2D  (  „  A 

*o3  yP’PJ) 


(7a) 


where 


CHEM2D 
ro 3 


(<p,p,t)  =  - 


d(P-L) 


dXc 


n-1 


(7b) 


is  the  ozone  photochemical  relaxation  time 
scale,  and  for  x0  we  relax  to  the  new  2D 
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Figure  8:  Ozone  photochemical  relaxation  times, 
in  days,  for  June  and  December  for  CHEM2D 
scheme  (top  row),  LINOZ  scheme  (middle  row) 
and  Cariolle-Deque  scheme  (bottom  row). 


ozone  climatology  discussed  in  section  3b. 
The  NOGAPS-ALPHA  lookup  tables  for  our 
CHEM2D  ozone  relaxation  times  are  plotted 
for  June  and  December  in  Figures  8a  and  8b, 
respectively.  The  time  scale  is  very  long  in  the 
lower  stratosphere  where  ozone  is  tracer-like, 
and  short  in  the  upper  stratosphere  where 
photochemistry  is  dominant.  The  rates  also  get 
very  long  at  all  heights  in  polar  night  where  no 
sunlight  is  available  to  drive  the 
photochemistry. 

LINOZ  scheme 

Two  linearized  ozone  photochemistry 
schemes  based  on  (5)-(6)  already  exist,  for 
which  we  have  obtained  the  4  governing 
linearized  photochemical  coefficients.  The 
first  is  the  so-called  LINOZ  scheme  of 
McLinden  et  al.  (2000).  The  LINOZ 
coefficients  were  computed  using  the 
University  of  California,  Irvine 
photochemistry  box  model.  The  scheme 
relaxes  via  (6)  to  reference  2D  climatologies 
for  ozone,  temperature  and  ozone  column 
abundances  published  in  technical  reports 
from  NASA’s  Airborne  Effects  of 
Stratospheric  Aircraft  (AESA),  which  we  also 
obtained.  The  LINOZ  ozone  photochemical 
time  scales  for  June  and  December  are  plotted 
in  Figures  8c  and  8d,  respectively,  and  span 
the  0.25-250  hPa  range  only:  above  and  below 
these  boundaries,  we  extrapolate  the  boundary 
values  upwards  and  downwards  respectively 
while  zeroing  out  the  other  extrapolated 
photochemical  coefficients,  yielding  a  simple 
relaxation  to  the  ozone  climatology. 

Cariolle-Deque  scheme 

A  full  version  of  (5)-(6)  was  developed  in 
the  mid-eighties  by  Cariolle  and  Deque  (1986) 
using  their  two-dimensional  chemistry  model, 
for  use  in  their  3D  climate  model.  That  same 
basic  scheme  has  recently  been  incorporated 
into  the  ECMWF  model,  with  linearized 
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Figure  9:  NOGAPS-ALPHA  ozone  mixing  ratio  forecasts  (in  ppmv)  over  Kiruna,  Sweden  for  runs 
initialized  on  1 1  January,  2003  at  0:00Z  and  progressing,  left  to  right  and  down,  in  +24  hour  forecast 
increments  from  analysis  to  +120  hour  (5-day)  forecasts.  Black  dashed  curves  shows  passive  O3,  and 
colored  curves  show  active  ozone  using  linearized  photochemistry  schemes  from  CHEM2D  (green), 
Cariolle-Deque  (blue),  and  LINOZ  (red).  Crosses  show  actual  ozone  measurements  on  these  days  from 
the  SAGE  III  satellite. 


photochemistry  coefficients  updated  using  the 
latest  version  of  that  2D  radiation-chemistry- 
dynamics  model  (Dethof  2003).  We  have 
obtained  those  latest  photochemical  lookup 
tables  and  implemented  them  into  a  third 
NOGAPS-ALPHA  ozone  photochemistry 
scheme  using  algorithms  similar  to  those  we 
developed  for  the  LINOZ  scheme. 

Figures  8e  and  8f  show  the  photochemical 
relaxation  times  for  June  and  December, 
respectively.  They  show  the  same  basic 
variations  as  the  CHEM2D  and  LINOZ 
coefficients.  The  photochemical  time  scales  in 
the  lower  stratosphere,  while  still  quite  long, 


are  somewhat  shorter  than  those  in  the  LINOZ 
and  CHEM2D  schemes. 

Intercomparison  Runs 

We  have  conducted  a  number  of 
intercomparison  runs,  both  single-column  tests 
and  full  NOGAPS-ALPHA  forecast  runs,  to 
assess  the  relative  merits  and  performance  of 
these  three  linearized  ozone  photochemistry 
schemes.  Some  of  this  work  will  be  reported 
more  fully  elsewhere  (McCormack  et  al. 
2004),  so  we  show  only  an  example  or  two 
here. 
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Figure  10:  ( P-L)0  coefficients  during  January  for  NOGAPS-ALPHA  linearized  03  photochemistry 
schemes  based  on:  (a)  CHEM2D  coefficients;  (b)  LINOZ  coefficients;  (c)  Cariolle-Deque  coefficients. 
CHEM2D  and  Cariolle-Deque  contours  are  in  units  of  10'2  ppmv  day'1,  while  LINOZ  values  are  in  ppmv 
day'1. 
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Figure  9  shows  a  series  of  NOGAPS- 
ALPHA  3D  ozone  profiles  over  Kiruna, 
Sweden  from  a  sequence  of  T79L54  forecast 
runs  initialized  on  11  February  2003  (0:00Z) 
using  MVOI/STRATOI  meteorological  initial 
conditions  and  3D  O3  initialization  fields  from 
the  GMAO.  At  this  time,  the  second  NASA 
SAGE  III  Ozone  Loss  and  Validation 
Experiment  (SOLVE  II)  was  underway,  with 
NASA’s  instrumented  DC-8  research  aircraft 
flying  out  of  Kiruna  (the  DC-8  flew  missions 
that  measured  ozone  and  related  parameters 
on  January  12,  14  and  16).  A  wealth  of  other 
data  were  acquired  during  SOLVE  II,  and 
many  global  models  were  run  for  this  period, 
and  so  this  SOLVE  II  period  focused  on 
Arctic  winter  represents  an  excellent  focus  for 
NOGAPS-ALPHA  runs  aimed  at  validating 
our  emerging  3D  prognostic  ozone  capability. 
On  the  other  hand,  many  of  these  high-latitude 
regions  during  SOLVE  II  are  in  the  polar 
night  and  so,  unless  this  air  was  transported  to 
lower  sunlit  latitudes,  they  may  not  undergo 
much  photochemistry,  making  relative 
evaluation  of  the  three  ozone  photochemistry 
schemes  a  little  more  difficult. 

Nonetheless,  Figure  9  reveals  some 
interesting  differences  among  predicted  ozone 
profiles  governed  by  each  of  the  three- 


different  linearized  photochemistry  schemes 
(colored  curves)  and  passively-advected  ozone 
(dashed  curve).  The  most  obvious  feature  is  a 
systematic  substantial  underprediction  of 
ozone  mixing  ratios  on  the  topside  of  the 
ozone  layer  using  the  LINOZ  scheme.  This  is 
a  recurrent  feature  of  nearly  all  of  our 
experiments  with  LINOZ,  both  in  3D  runs 
using  NOGAPS-ALPHA  and  in  simpler 
offline  single-column  tests.  McLinden  et  al. 
(2000)  comment  on  this  propensity,  attributing 
it  to  the  well  known  “ozone  deficit”  problem 
encountered  in  upper  stratospheric  ozone 
modeling:  however,  model  ozone  deficits  are 
generally  smaller  than  this,  and  newer  data 
and  models  now  seem  to  have  largely 
eliminated  this  discrepancy  between 
observations  and  models  (Siskind  et  al.  1998; 
Groos  et  al.  1999). 

Instead,  these  underestimates  seem  to 
result  primarily  from  excessively  large 
negative  equilibrium  production/loss 

coefficients,  (P-L)o,  in  the  LINOZ  scheme  in 
the  upper  stratosphere.  Figure  10  compares 
(P-L)o  values  among  the  three  schemes  for 
January.  The  CHEM2D  and  Cariolle-Deque 
values  are  comparable  in  magnitude  and  have 
similar  latitude -height  distributions. 

Furthermore,  the  NCEP  Global  Forecasting 
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System  (GFS)  uses  (P-L)o  values  from  the 
Goddard  two-dimensional  model  (Fleming  et 
al.  2002)  as  their  current  ozone 
photochemistry  scheme,  and  those  values  too 
look  very  similar  to  those  in  Figures  10a  and 
10c.  The  LINOZ  values,  on  the  other  hand, 
look  very  different,  with  loss  at  all  latitudes  in 
the  upper  stratosphere  and  absolute  values  ~1- 
2  orders  of  magnitude  larger  than  the  other 
two  schemes.  We  speculate  that  this  may 
possibly  have  resulted  in  part  from  use  of 
ozone  rather  than  scaled  odd-oxygen  rates  in 
the  calculations  here. 

Whatever  the  reasons,  these  discrepancies 
make  the  current  LINOZ  scheme 
unsatisfactory  as  an  ozone  photochemistry 
parameterization  for  NOGAPS-ALPHA.  For 
example,  since  SW  O3  heating  peaks  at  ~1  hPa 
(Figure  2),  underpredictions  of  O3  mixing 
ratios  of  the  magnitude  noted  for  LINOZ  in 
Figure  9  would  lead  to  massively 
underestimated  stratospheric  heating  were  we 
to  use  3D  LINOZ  O3  fields  in  our  radiation 
calculations. 

Figure  9  also  plots  measured  ozone 
profiles  over  this  location  for  each  forecast 
day  from  the  SAGE  III  satellite  instrument 
(Poole  et  al.  2003).  While  LINOZ  clearly 
underpredicts  these  data,  the  other  two 
schemes  (CHEM2D,  Cariolle-Deque)  appear 
to  reproduce  the  SAGE  III  observations  fairly 
well,  and  in  this  case  it  is  difficult  to 
distinguish  any  differences  in  skill  between 
the  two.  Indeed,  even  passive  ozone  is  a 
reasonable  fit  to  these  data,  which  is  not 
surprising  given  the  long  photochemical  time 
scales  to  be  expected  in  polar  night  (Figure  8). 
There  is  some  suggestion  that  ozone  in  the 
lower  stratosphere  (50-100  hPa)  may  be 
underpredicted:  given  the  long  photochemical 
time  scales  at  all  latitudes  at  50-100  hPa  (see 
Figure  8),  this  difference  may  reflect  a  slight 
error  in  our  ozone  initialization  fields  for  this 
period.  McCormack  et  al.  (2004)  provide  a 
thorough  study  of  NOGAPS-ALPHA  ozone 


forecasts  for  various  periods  during  SOLVE 
II. 

Another  useful  NWP  test  is  to  compute 
synoptic  maps  of  forecast  total  ozone,  which 
is  important,  for  example,  in  forecasting 
surface  UV  levels  (Long  2003).  Figure  11 
shows  sample  NOGAPS-ALPHA  total  ozone 
forecasts  initialized  on  20  September  2002 
during  the  unprecedented  splitting  of  the 
Antarctic  Ozone  Hole  due  to  a  major 
stratospheric  warming  event  (see,  e.g.,  Allen 
et  al.  2003,  2004).  We  see  that  NOGAPS- 
ALPHA  captures  the  vortex  splitting  event 
and  the  separation  of  the  Ozone  Hole  into  two 
separate  pieces,  with  the  general  morphology 
and  dynamic  range  of  the  total  ozone  values 
quite  similar  to  measurements  from  the  Total 
Ozone  Mapping  Spectrometer  (TOMS). 

5.  OUTLOOK 

As  some  of  the  features  of  NOGAPS- 
ALPHA  that  we  have  introduced  in  this  paper 
evolve  and  mature,  they  will  be  progressively 
transitioned  to  operations  at  FNMOC.  The 
new  high-altitude  prognostic  capabilities  they 
will  provide  to  NOGAPS  should  reap  benefits 
in  overall  NWP  skill  and  provide  new  forecast 
products  of  direct  relevance  to  the  Navy  and 
the  Department  of  Defense  (DoD)  generally. 
We  conclude  by  discussing  some  of  these 
potential  benefits  and  new  uses. 

Efforts  are  already  underway  to  assimilate 
satellite  radiances  directly  using  NAVDAS. 
The  high-altitude  capabilities  of  NOGAPS- 
ALPHA  can  yield  considerable  improvements 
in  these  assimilations.  Prognostic  3D  ozone, 
for  example,  can  help  correct  complex  biases 
in  some  LW  radiance  channels  due  to  ozone 
absorption  (Derber  and  Wu  1998).  Improved 
operational  assimilation  should  then  lead  to 
improvements  in  NWP  skill  at  all  altitudes. 

There  is  an  emerging  realization  that 
important  regional  winter  weather  influences, 
such  as  the  North  Atlantic  Oscillation  (NAO), 
are  (in  part)  manifestations  of  a  global 
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Figure  11:  Top  row:  total  ozone  (in  Dobson  units)  over  the  South  Pole  derived  from  NOGAPS-ALPHA 
T79L54  forecasts  initialized  on  20  September  2003  (12:00  Z)  using  3D  ozone  initialized  from  GMAO 
analyses  and  updated  using  the  Cariolle-Deque  photochemistry  scheme.  Results  are  plotted  at  forecast 
day  1  (21  September),  day  3  (23  September)  and  day  5  (25  September).  Beneath  each  plot  is  the 
corresponding  total  ozone  measured  on  this  day  by  TOMS. 


atmospheric  phenomenon  known  assortedly  as 
the  Arctic  Oscillation  (AO)  or  the 
Northem/Southem  Annular  Mode 
(NAM/SAM).  Recent  research  has  shown  that 
anomalies  in  the  strength  of  the  AO  often  form 
in  the  upper  stratosphere  and  then  propagate 
downwards  over  a  series  of  days,  sometimes  to 
the  surface  to  change  the  state  of  the  surface 
AO,  NAO  and  hence  global  weather  patterns 
(Baldwin  and  Dunkerton  2001;  Thompson  and 
Wallace  2002;  Baldwin  et  al.  2003;  Charlton 
et  al.  2003).  Only  a  high-altitude  global  NWP 
model  like  NOGAPS-ALPHA  is  capable  of 
capturing  these  kinds  of  global  AO/NAO 
dynamics  and  their  effect  on  NWP. 

Entirely  new  high-altitude  products  will  be 
provided  by  a  NOGAPS-ALPHA  prototype 
running  operationally,  and  these  should  serve 


many  areas  of  the  military,  which  currently 
lack  adequate  forecast  guidance  at 
stratospheric  and  mesospheric  altitudes. 
Global  stratospheric  forecasts  from  NOGAPS- 
ALPHA  would  support  the  rapidly  expanding 
fleet  of  high-altitude  long-endurance  (HALE) 
UAVs,  such  as  Global  Hawk,  which  are 
expensive  airframes  that  cruise  in  the  middle 
atmosphere,  are  structurally/aerodynamically 
vulnerable  to  atmospheric  loads  (Eckermann  et 
al.  2004),  and  have  almost  global  range.  High- 
altitude  forecasts  are  also  needed  to  assess  (in 
real-time)  loads  for  boost-phase  missile 
defense  (Kleppner  et  al.  2004)  and  refractivity/ 
scintillation  forecasts  for  the  Airborne  Laser 
(ABL)  (Hecht  2004).  Similar  issues  are 
relevant  for  DoD  spacecraft  launch  and  reentry 
projections,  and  real-time  signal  correction 
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between  DoD  satellites  and  forces/targets  on 
the  ground  or  in  the  air  (e.g.,  UAVs).  High- 
yield  nuclear  detonations  inject  significant 
volumes  of  radioactive  debris  well  into  the 
stratosphere,  where  it  can  be  rapidly  dispersed 
and  transported  globally  (see,  e.g.,  Figure  6). 
The  prognostic  chemical  capabilities 
developed  here  for  NOGAPS-ALPHA  can  in 
principle  be  extended  to  any  injected  tracer, 
and  thus  could  be  applied  at  short  notice  to 
provide  inline  global  forecasts  of  projected 
dispersal  and  fallout  of  such  hazardous 
material  (e.g.,  Feely  et  al.  1966a  1966b; 
Rehfeld  and  Heimann  1995).  This  internal 
chemical  transport  functionality  in  NOGAPS- 
ALPHA  is  particularly  attractive  for  this  task, 
since  recent  work  has  shown  that  transport 
using  the  internal  model  dynamics  of  a  GSM 
is  superior  compared  that  from  offline 
transport  models  using  output  from 
assimilation  systems  (e.g.,  Schoeberl  et  al. 
2003). 
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